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corresponding  more  direct  method  to  compute  the  natural  frequencies  in  exterior 
scattering  problems.* 

\ 

To  partially  answer  this  question  a  geometric  ray  optics  method,  which  because 
of  its  asymptotic  nature  is  particularly  suited  to  compute  the  higher  order 
resonances ,  is  described. 


The  idea  is,  in  general  terms,  to  consider  a  smooth,  convex  object  with  a  sur¬ 
face  impedance  boundary  condition.  Complex  resonant  frequencies  are  computed 
from  closed  path  integrals  over  the  surface.  It  is  hypothesized  that  the 
paths  Ti  represent  surface  geodesics  whose  definition  includes  electromagnetic 
inertial  effects. 


The  electrical  path  length  depends  upon  the  lochl  principal  radii  of  curvature 
of  the  surface.  The  primary  result  of  this  work  is  that  a  generalized  WKB 
method  can  be  employed  to  account  for  both  the  curvature  and  Impedance  boundary 
conditions. 
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ABSTRACT 


Th*  representation  of  the  transient  electromagnetic  response 
of  finite  size,  smooth,  perfectly  conducting  objects  in  terms  of  a 
complex  exponential  series  is  the  central  ingredient  of  the  singular¬ 
ity  expansion  method  (SEM) .  The  exponential  terms  correspond  to  che 
complex  natural  frequencies  associated  with  the  object  geometry. 

That  such  a  simple  series  can  predict  the  force  free  response  of 
complicated  objects  begs  the  question  "is  there  not  a  corresponding 
more  direct  method  to  compute  the  natural  frequencies  in  exterior 
scattering  problems." 

To  partially  answer  this  question  a  geometric  ray  optics 
method,  which  because  of  its  asymptotic  nature  is  particularly  suited 
to  compute  the  higher  order  resonances,  is  described. 

The  idea  is,  in  general  terms,  to  consider  a  smooth,  convex 
object  with  a  surface  impedance  boundary  condition.  Complex  resonant 
frequencies  are  computed  from  closed  path  integrals  over  the  surface. 
It  is  hypothesized  that  the  paths  represent  surface  geodesics 
whose  definition  includes  electromagnetic  inertial  effects. 

The  electrical  path  length  depends  upon  the  local  principal 
radii  of  curvature  of  the  surface.  The  primary  result  of  this  work 
is  that  a  generalized  WKB  method  can  be  employed  to  account  for  both 
the  curvature  and  impedance  boundary  conditions. 


The  mechod  when  applied  to  a  sphere  is  shown  to  reduce  to 
the  well  known  uniform  asymptotic  expansion  of  the  spherical  Hankel 
functions.  In  this  example  a  comparison  of  the  asymptotic  and 
exact  results  for  the  natural  oscillation  frequencies  of  a  sphere 
is  given. 


v 


Introduction 


The  complex  poles  in  Che  s  domain  which  determine  Che  complex 
resonant  frequencies  of  finite  convex  scattering  geometries  have  been 
determined  numerically  [1].  From  such  studies  it  has  been  observed 
that  the  complex  frequencies  or  pole  locations  depend  only  upon  the 
scattering  obstacle  and  not  upon  the  incident  waveform.  In  ref¬ 
erence  [1]  ,  Tesche  observed  that  the  approximate  distance  between 
adjacent  poles  in  the  s  plane  have  the  property  that 

Im  (A  s)  *  tt  c/L  (1) 

where  L  is  a  characteristic  length  of  the  body  and  c  is  the  speed  of 
light.  Physically,  one  expects  that  the  exterior  frequencies  must 
be  complex  to  account  for  ray  divergence  and  bending  losses. 

These  observations  suggest  that  the  exterior  resonance  phenom¬ 
enon  is  a  form  of  damped  periodic  motion.  Similar  systems  are 
cavity  resonance  and  Bohr  orbits.  Thus  (1)  can  be  considered  to  be 
related  to  a  phase  reinforcement  condition  which  requires  that  closed 
ray  paths  around  the  object  corresponding  to  resonant  frequencies 
effectively  contain  an  integer  number  of  wavelengths.  This,  of 
course,  is  the  principle  of  Bohr  orbits  of  the  Old  Quantum  Theory  [2]. 

An  asymptotic  geometric  theory  should  apply  to  electrically 
large  objects  and  improve  as  the  object  dimensions,  scaled  in  wave¬ 
lengths,  or  mode  index  n  becomes  large.  The  quasi-periodic  nature  of 
observed  transient  responses  from  such  objects  suggest  the  candidate 
theories  of  transverse  resonances  [3],  the  Wentzel-Kramers-Brillouin  or 
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WKB  approximation  [4],  and  a  suitable  modification  of  the  Bohr- 
Sommerfeld-Wilson  condition  of  the  Old  Quantum  Theory  [2]. 

Theory 

To  motivate  the  development  to  follow,  it  is  useful  to  out¬ 
line  the  elegant  geometric  treatment  of  Interior  resonances  of 
Keller  and  Rubinow  [5].  Many  of  their  basic  ideas  and  topological 
concepts,  as  well  as  their  comparative  accuracy,  can  be  expected  to 
carry  over  to  the  more  difficult  external  problem.  Their  idea  is 
to  solve  the  scalar  wave  equation 

(72  +  k2n2)  <Kx)  =  0  (2) 

on  the  interior  of  a  smooth  closed  surface  S  where  is  assumed  to 
satisfy  Neumann  or  Dirichlet  conditions  on  S.  If  we  assume  a 
solution  of  the  fora 

<Kx)  »  A  (x)eikoS(x)  (3) 

it  follows  from  substituting  (3)  into  (2)  that  the  phase  function 
S(x)  and  the  amplitude  function  A(x)  satisfy  the  eikonal  and  trans¬ 
port  equations 

(V  S)2  -  n2  (4) 

V2  S  +  2  7  S  •  In  A  -  0  .  (5) 

If  the  local  ray  trajectory  has  a  unit  tangent  vector  £  then  the  local 

wave  vector  k  is  given  by 

k  -  kQ  7  S  ,  7  S  •  £  n(x)  (6) 
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The  resonance  condition  is 


dl 


2tt  (n  +  m/4  +  b/2) 


(7) 


where  n,  m,  and  b  are  integers.  In  particular,  m  and  b  are  the  number 
of  times  the  ray  hits  the  caustic  or  boundary  respectively.  (This  is 
for  Dirichlet  conditions  on  S;  for  the  Neumann  condition,  b  is 
identically  zero.)  Condition  (7)  holds  for  every  independent  closed 
curve  on  the  covering  space.  Keller  and  Rubinow  consider  the  example 
of  a  resonant  ray  inside  a  perfectly  conducting  circular  surface  of 
radius  a.  Using  geometric  optics,  one  finds  that  there  exists  a 
caustic  surface  of  radius  a^  <  a  as  shown  in  Figure  1. 


Caustic 
(radius  aQ) 

Ray 

Boundary 
( radius  a) 


Figure  1.  Interior  Resonant  Ray 
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If  one  accounts  for  the  number  of  ray  congruencies,  and  the 
independent  closed  curves,  by  successive  elimination  of  the  parameters 
such  as  the  caustic  radius  aQ  in  Figure  1,  the  desired  wave  number  k 
and  hence  the  complex  frequency  3  can  be  obtained. 

The  independent  paths  on  the  covering  cannot  be  deformed  into 
each  other  without  crossing  a  singularity  of  the  field.  Thus, 

Equation  (7)  is  actually  a  residue  theorem  in  disguise.  This  charac¬ 
terization  of  the  resonance  condition  can  be  used  directly  in  target 
discrimination  as  we  will  discuss  in  conclusion  of  the  paper. 

In  the  example  of  Figure  1  there  are  two  independent  paths 
in  the  covering  labeled  1  and  2  in  Figure  2. 


Caustic 
(radius  a0) 


k  2ira  ■  27rm  (8) 
oo 


Ray  -  -  Boundary 
(radius  a) 


k  (2  /z2-a2  -  2a  cos_1(A))  (9) 

0  0  0 

*  2'!T(n  +  3/4) 


Figure  2.  The  Independent  Paths 


Eliminate  aQ  from  Equations  (8)  and  9)  and  solve  for  ka. 

The  resulting  expression  for  ka  is  an  asymptotic  formula  for  j  , 

*  *  nm 

the  nth  zero  of  the  Bessel  function  J  . 

m 

The  following  table  from  Reference  2  compares  approximate 
and  exact  zeros. 


Table  1.  Interior  Resonance  Comparison 


m 

n 

Approximate 

Exact 

Fractional  Error 

0 

1 

2.356 

2.405 

0.0204 

0 

2 

5.498 

5.520 

0.0040 

Q 

3 

8.639 

3.654 

0.0017 

0 

4 

11.781 

11.792 

0.0009 

1 

1 

3.795 

3.802 

0.0097 

1 

2 

6.997 

7.016 

0.0027 

1 

3 

10.161 

10.173 

0.0012 

1 

• 

4 

13.311 

• 

13.324 

0.0010 

• 

• 

• 

• 

• 

• 

• 

• 

As  can  be  seen,  the  asymptotic  solution  is  surprisingly  accurate 
even  for  the  lowest  order  modes. 

In  the  exterior  problem  the  situation  becomes  more  complicated. 
The  role  of  the  caustic  and  the  boundary  are  interchanged  and  the  ray 
paths  become  curved.  In  Figure  3,  evanescent  rays  are  shown  trapped 
between  the  surface  Sq  and  their  respective  caustics.  When  curvature  is 
encountered,  some  of  the  energy  radiates  away  from  the  object. 
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Caustic 

Surface 


Figure  3.  Exterior  Resonance  Ray  Orbits 

It  is  easy  to  appreciate  that  the  corresponding  eigenvalues 

in  the  exterior  problem  are  complex.  The  imaginary  part  of  the  eigen- 
frequencies  sn  associated  with  scalar  transient  f(t)  are  associated 

with  the  object  circumference.  Let  f(t)  be  written  as  the  inverse 

Laplace  transform  of  F(s). 

rc+i® 

f(t)  *  27Ti  F(s)  eSt  ds 

J  c-i® 

Let  s  be  a  simple  pole  of  F(s)  in  the  left  half  plane  of  the  s 
n 

domain ,  i . e . , 

s  ■  o  +  i  co  ,  a  <  0. 
n  n  n  n 
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The  real  part  0n  is  physically  associated  with  radiation  damping  and 
ray  tube  spreading  upon  reflection  from  the  convex  surface.  The  Idea 
Is  demonstrated  In  Figure  4. 


Figure  4.  Heuristic  Definition  of  Caustic 

The  caustic  surface  for  the  external  problem  can  heuristically 
be  defined  as  the  locus  where  the  local  wave  front  tangential 
velocity  exceeds  that  vacuum  velocity  c.  The  energy  must,  therefore, 
detach  from  the  packet.  This  phenomenon  is  associated  with  the  Inertia 
of  the  electromagnetic  field:  there  is  a  cause  and  effect  relation¬ 
ship  between  local  spatial  curvature  and  electromagnetic  energy 
density. 

Such  a  relationship  was  obtained  by  Einstein  in  his  theory 
of  General  Relativity  in  the  form  of  a  conservation  law  involving  the 
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covariant  derivative  of  the  electromagnetic  energy  momentum  tensor 
Ty\>  The  mathematical  statement  is  given  by 

3/3x  (/£  TV  x)  -  \  (-%)  /g  Tuv  -  0  (10) 

V  ^  2  3xX 

where  gUV  is  the  space  time  metric  tensor  and  g  ■  det  (glJV)  [6,7]. 

It  has  been  shown  by  Choudhary  and  Felsen  [8]  that  geometric 
ray  tracing  in  evanescent  regions  is  further  complicated  by  the  non¬ 
congruence  of  the  phase  propagation  paths  of  simple  geometric  theory 

and  the  power  flow  trajectories. 

To  obtain  a  geometric  theory  for  the  exterior  eigenfrequencies, 

we  postulate  that  the  periodic  damped  motion  be  described  by  a 
suitable  modification  of  the  Bohr-Sommerfeld  quantization  rule  for  the 

action  integral  of  the  "Old  Quantum  Theory,"  [2,9,10], 

j>  p^  dqA  *  n^h  ,  i  *  1,2,3...,  n^  *  1,2,3.  (11) 

Here,  p^,  q^  are  the  canonically  conjugate  momenta  and  coordinates  as 
defined  by  the  equations  of  motion  in  the  Hamiltonian  formalism.  An 
important  point  is  that  for  each  independent  momentum  coordinate  set 
(Pf,  qi)  ,  i  *  1,2,...  there  is  a  corresponding  quantization  or 
resonance  condition.  In  the  classical  physics  problem  of  determining 
the  resonant  frequencies  of  an  electromagnetic  cavity  oscillator,  Keller 
and  Rubinow  [3]  showed  that  Equation  (7)  can  be  used  to  compute  the 
associated  eigenvalue  spectrum. 

We  want  to  show  that  analogous  independent  closed  curves 

r 

p  . 


p  -  1,2,3. . . 


can  also  be  determined  and  that  the  independent 


momentum  -  coordinate  defined  action  integrals  for  the  exterior  problem 
determine  the  eigenvalue  spectrum. 

The  Bohr-Sommerfeld  rule  was  developed  for  the  determination 
of  electron  orbits  and  associated  energy  eigenvalues  of  atomic  physics. 
However,  the  wave  nature  of  particles  is  manifest  according  to 
de  Broglie's  relation  that  a  particle  with  momentum  p  has  associated 
wave  properties  with  wavelength  X  [7].  This  relation  is 

X  »  h/p  or  k.  *  Q  »  2irp  (12) 

A  h 

Thus,  the  wavenumber  momentum  relation  corresponding  to  (1)  is 

j  k±  dqt  *  2im1  i  «  1,2,3...  (13) 

r±  n±  ■  0,1,2,3. . . 

We  postulate  that  (13)  with  proper  interpretation  determines  asymptot¬ 
ically  the  complex  resonant  frequencies  of  an  exterior  scatterer.  In 
the  original  development  of  the  Old  Quantum  Theory,  the  correspondence 
principle  was  used  to  bridge  the  gap  between  classical  and  quantum 
mechanical  domains.  Then,  the  classical  limit  is  approached  as  the 
integer  n^  in  Equation  (11)  becomes  large.  Similarly,  the  asymptotic 
optical  limit  for  electromagnetic  analogue  in  Equation  (13)  is 
approached  when  the  corresponding  n^  becomes  large.  To  apply  these 
concepts  we  first  develop  the  relevant  scattering  problem. 

Our  boundary  value  problem  is  as  follows.  We  wish  to  solve 
the  scalar  wave  Equation  (2)  for  ^(x)  in  the  region  exterior  to  the 


surface  SQ  as  shown  in  Figure  3.  The  boundary  conditions  on  ’P(x) 
is  an  impedance  condition  of  the  form 

-  0  (14) 

X£  S 

0 

where  q  ■  -i  ^  s  z  A  ,  z  -  120  tr  ft,  and 
ooo 

A  *  normalized  surface  impedance  at  grazing  incidence 

The  same  type  of  asymptotic  solution  described  in  Equations  (3) 
through  (7)  are  again  used.  The  equation  corresponding  to  (13)  for  the 
exterior  ray  orbits  over  an  impedance  surface  is 

kQ  |  V  S  *  dl  -  2ir  (nt  -  -f  -  *4)  (15) 

'r±  i  *  1,2,3  ... 

The  interpretation  of  Equation  (15)  is  the  same  as  (13)  where  now  the 
phase  ^  is  given  by 

♦i*2?  “s  <n.  r-|ji  (i6) 

where  A  is  the  normalized  outward  pointing  surface  impedance  intro¬ 
duced  in  Equation  (14).  Thus 

arg  (D  ■  Im(ln  (?)) 

and  the  paths  F^  are  yet  to  be  determined.  The  caustic  in  the  exterior 
problem,  as  shown  in  Figure  3,  separates  the  dark  from  the  light  side. 
The  dark  side  is  closer  to  the  surface  S  ;  outside  the  caustic  the 


<§  ♦  «» 
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energy  radiates  away  from  the  surface.  The  caustic  or  turning  point 
then  separates  the  evanescent  and  radiating  regions.  The  fact  that 
the  phase  changes  it/ 2  radians  at  the  caustic  can  be  seen  from 
Equation  (16)  if  the  appropriate  normalized  impedance  A  ■  i  is  used. 

As  described  by  Howard  [lXI .  this  caustic  phase  advance  is 
central  to  the  bending  loss  associated  with  open  waveguides.  In 
addition  to  the  bending  losses  and  the  ohmic  losses  in  the  scatter¬ 
ing  object,  there  is  a  third  loss  mechanism  in  exterior  problems. 

This  is  the  ray  spreading  caused  by  reflection  off  the  convex 
surface  Sq.  One  way  to  account  for  this  ray  tube  divergence  is  to 
place  the  amplitude  function  A(x)  into  the  exponent  and  then  replace 
Vs  in  Equation  (15)  by 

VS  -  i/k  V  In  A  (17) 

o 

The  dependence  of  the  amplitude  function  A  near  a  surface  such  as  Sq 
has  been  carried  out  by  Kouyoumj ian  and  Pathak  [ 12 , 13 ] . 

APPLICATION  TO  A  SPHERE 

The  exterior  resonances  of  a  sphere  have  been  thoroughly  inves¬ 
tigated  beginning  with  Thompson's  1884  treatment  [14].  A  more  modem 
discussion  can  be  found  in  Stratton's  book  [15].  Numerical  results  have 
been  presented  by  Martinez  et  al  [16].  In  these  references  it  is  shown 
that  the  complex  resonances  for  the  sphere  are  associated  with  the  com¬ 
plex  zeros  Pnm  of  the  spherical  outward  radiating  Hankel  function 
h^^ .  In  the  case  of  perfect  conductivity  the  mode  equations  are 

H  Type  Modes  (Horizontal  Polarization,  q  ■*  *) 

hf1)(P  )  -  0  (18) 

n  nm 


ii 


(19) 


E  Type  Modes  (Vertical  Polarization,  q  ♦  0) . 

(ph<1}(p))'  -  0 

The  poles  s  in  the  Lapace  s  domain  which  corresponds  to  the  complex 
nm 

resonant  frequencies  are  given  by 

P„m 

s  -  -i  —  ,  m  -  1,2,3  ...  (20) 

a^F 

where  a  is  the  sphere  radius. 

Because  the  spherical  Hankel  function  of  order  n  is  a  polynomial  of 
order  n  multiplied  by  an  exponential  function,  the  natural  frequencies 
for  a  sphere  are  particularly  easy  to  compute. 

The  geometric  theory  as  represented  by  the  generalized  Bohr- 
Sommerfeld-Wilson  condition  (15)  should  asymptotically  produce  the 
complex  frequencies  as  defined  by  Equations  (18),  (19)  and  (20). 

For  the  general  case  of  an  arbitrarily  shaped  convex  scattering 
body,  a  major  difficulty  in  applying  Equation  (15)  is  in  the  determination 
of  the  fundamental  paths  r  .  It  is  anticipated  that  the  required  geo¬ 
desic  paths  will  be  best  computed  using  local  geodesic  coordinates  on 
the  surface  of  SQ  as  defined,  for  example  by  Struik  [17]  or  O'Neill  [18]. 
These  local  solutions  then  are  integrated  over  the  global  geodesic  paths 
to  determine  the  resonant  frequencies. 

For  a  sphere,  the  local  solution  can  be  trivially  translated 
into  the  global  solution.  A  coordinate  independent  general  way  to 
approach  the  problem  is  to  compute  the  ray  divergence  contribution  to 
the  loss  term  in  Equation  (17)  through  the  geometrical  formulation  of 
reflection  from  curved  surfaces  as  developed  for  example  by 
Kouyoumj ian  and  Pathak  [ 12 ] . 
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Keller  has  carried  out  this  approach  for  asymptotic  solutions 
to  the  Schrodinger  equation  [2].  He  found  that  the  amplitude  contribu¬ 
tion  to  the  phase  (as  determined  by  Equation  (17)  here,  for  example) 

t 

gives  rise  to  half-integer  quantum  numbers. 

At  this  stage  of  progress  in  our  work,  a  heuristic  approach 
to  the  geometric  computation  of  the  exterior  resonances  of  a  metallic 
sphere  is  given. 

The  convex  scattering  surface,  as  shown  in  Figure  3,  can  be 
defined  mathematically  by  an  equation  of  the  form 

f(x,y,z)  *  0  (21) 

In  a  neighborhood  of  the  point  (0,0,0)  it  is  possible  to  represent  f 
in  the  form 

f  »  z  +  1/2  (ax2  +  28xy  +  Yy2)  (22) 


where  the  convexity  of  the  surface  and  fixing  the  orientation  of 
positive  z  to  the  convex  side  of  f  requires  that 


a  >  0,  8  >  0,  ay  -  S2  >  0 


(23) 


The  coefficients  a,  8,  y  are  given  by 

32z 


a  «  - 


2 

a  z 


ctx 


*  -  y  -  0*  S’"3x3y 


x  -  y  -  0 


0  (24) 


so  that  representation  (22)  is  just  a  McLaurln  series  through  second 
order  terms  about  the  point  x  *  y  ■  0  of  the  surface  (21).  At  the 
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point  (0,0,0)  on  the  surface  the  normal  curvature  in  the  x  direction 
is  a  ■  1/a  and  the  analogous  y  quantity  is  b  »  1/8. 

To  model  a  convex  body  with  two  principal  radii  of  curvature 

we  choose  the  system  defined  by  the  arc  length  formula 

2  2  2  2  2  2 
ds  -  p  d0  +  b  d<T  +  dp  ‘ 


The  Laplacian  operator  is  determined  to  be 


T7 a,.,  1  f  3  ,b  ,  3  ,p  ,  3  .  .  , 

7  *  ‘  Sb  >39  +  1?  <b  W  *  Ip  <Pb 


So' 


(2i 


Consider  ray  solutions  propagating  in  the  direction  of  Increasing  6. 
We  transform  out  the  transverse  $  dependence.  Thus,  let 


’P(p,M) 


r 

*  - JX 


'P  (p,0,t)  elk<(>t  dt 


(2f 


The  wave  equation 


(V2  +  k2)  'p 


becomes  upon  using  (25)  and  (26), 


F3p  (p  +  p^I^+  k*(1  '  c 


Assume  a  product  solution  of  the  form 


7b2)  it  •  0 


iykaQ 


i>  (p,9,t)  -  R(p)  e 


where  a  is  the  local  radius  of  curvature  in  the  ray  direction,  and  y  is 
a  normalized  propagation  constant  to  be  determined.  The  ordinary  dif¬ 
ferential  equation  for  R(p)  is  thus 

P  dp  (P  dP}  +  <  (kP>2  -  kV  ($2-  Y2(ka)2)R(p)-  0  (29) 

Let  x  -  kp(l  -  (r))1^2  »  v  *  Yka  ,  u(x)  -  R  (p) 

0 

Then,  the  equation  for  u(x)  is 


x  (x  u’OOr  +  (x2  -  \>2)  u(x)  -  0 


In  equation  (30),  make  the  change  of  variable 


y  -  In  (x/x^)  ,  xa  ■  x  |p  ,  a  *  kaa  , 


f(y)  *  u(x) 


a  -  (1  -  (t/b)z) 


2,1/2 


The  resulting  differential  equation  for  f(y)  is 
— ' M  +  ft2n2  (y)  f  (y)  -  0 


where 


n2(y)  -  Cl2  e2?  “  Y2 


,  J2  *  ka 


The  radial  Bohr-Sommerfeld  resonance  condition  for  equation  (31) ,  which 
corresponds  to  a  ray  trip  up  and  back  between  sq  and  the  caustic  in 
Figure  3,  becomes 

2  8  fy°  n(y)dy  -  -  it/2  (32) 


♦ 


Here  we  have  chosen  ■  1  in  Equation  (15)  and  yQ  is  the  turning 
point  of  Equation  (31),  i.e.,  n(yQ)  *  0.  Also  in  Equation  (32), 
as  defined  by  Equation  (16)  is  zero.  This  is  the  case  for 
vertical  polarized  E  field  on  a  perfectly  conducting  sphere.  In 
Equation  (32),  we  make  the  change  of  variable 


T  *  e  <yo  ^  and  define  z  *  —  . 

y 

Then  Equation  (32)  becomes 

2ifrr  (  dx  -  -IT/2  . 


(33) 


Olver  has  made  use  of  this  Integral  extensively  in  his  work  on  special 
functions  [19,20].  He  defines  the  implicit  relationship  z(C)  as 

rz  7  .2%1/2 


2/3  ; 


3/2 


(- 


1  -  tl) 


dt 


(34) 


In  our  application  X,  is  given  and  z  is  to  be  computed.  This  is  a  complex 
valued  transcendental  equation  in  which  it  is  quite  useful  to  introduce 
an  intermediate  change  of  variable 

z  -  secho  (35) 

since  then  it  can  be  shown  that 

2/3  ;3/,2»  a  -  tanh3  (36) 

Details  can  be  found  in  Appendix  A. 

Equation  (33)  then  is  the  radial  resonance  condition.  For  the 
sphere,  the  angular  resonance  condition  from  Equation  (15)  and  our 
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assumed  angular  dependence  as  given  in  relation  (28)  is 


2iryft  -  2-rn  n  -  1,2,3,  ft  -  kQa  (37) 

since  the  ray  orbits  are  obviously  great  circles.  (All  great  circles 
on  a  sphere  are  geodesics.)  Combining  Equations  (33),  (34)  and  (37) 
yields  the  mode  equation 

— g  -  -ff/2  (mod  2tt)  .  (38) 

This  is  a  phase  requirement  so  that  £(z)  as  defined  by  Equation  (38) 
is  multivalued.  Thus,  from  Equation  (38)  it  follows 


2/3 


r  3/2 


-  g  (4m  -  1) 


m  ■  1,2,3. . . 


(39) 


On  the  other  hand,  the  uniform  asymptotic  formula  for  the  roots  of  the 
spherical  Hankel  function  h^  (z)  are  given  by  [2l] 


n 2/3  ?  (Z  )  ei27T/3  “  am  (40) 

nm 

where  a^  is  the  m^  root  of  the  Airy  function  Ai  [21].  The  zeros  of 
Ai  are  along  the  negative  real  axis.  The  lead  term  of  the  asymptotic 
expansion  for  the  zeros  of  Ai  is  (see  [21],  p  450) 


a 

m 


-  -  ( 3tt  (4m  -  l)/8)2/3 


m  »  1,2,3  . . . 


(41) 


If  expression  (41)  is  substituted  into  (40),  the  result  is  identical  to 
the  mode  equation  (39).  Thus,  ve  have  shown  that  the  modified  Bohr- 
Sommerfeld-Wilson  quantum  conditions  when  applied  to  the  exterior 
resonances  of  a  sphere  result  in  a  mode  equation  which  is  identical  to 
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the  uniform  asymptotic  formulas  for  the  complex  zeros  of  the  Hankel 
functions  as  given  by  Abramowitz  and  Stegun  [21]  or  Olver  [19]. 

To  give  an  indication  of  the  accuracy  of  the  method,  a  table 
of  the  exact  (Z^)  and  asymptotic  (2^  as  computed  using  Equation  (39) 
and  results  of  Appendix  A)  complex  zeros  of  h^^ (z)  is  given. 

Table  2.  Exterior  Resonance  Comparison  of  Exact  (Z  .  )  and 

nic 


Asymptotic  (Znk)  Zeros  of  h4^  (z). 


k 

(n  =  4) 

Znk 

Znk 

1 

2.6574, 

-2.1038 

2.6539, 

-2.1035 

2 

0.8672, 

-2.8962 

0.8628, 

-2.8950 

3 

-0.8672, 

-2.8962 

-0.8741, 

-2.8927 

4 

-2.6574, 

-2.1038 

-2.6740, 

-2.0892 

Again,  as  in  the  interior  resonance  comparison  of  Table  1, 
the  geometric  results  are  surprisingly  accurate. 

Discussions  and  Recommendations 

In  conclusion,  we  have  shown  that  a  modification  of  the  Bohr- 
Sommerfeld-Wilson  quantum  condition  can  predict  exterior  resonance 
frequencies.  The  method  is  geometric  and  hence  applies  to  non- 
separable  geometries.  It  has  the  potential  to  order  the  pole  trajec¬ 
tory  plots  of  SEM  into  radial  and  angular  families. 

It  is  important  to  realize  that  the  poles  we  have  determined 
do  not  correspond  to  creeping  waves.  Thus,  if  s  ^  *s  t*ie  position 
of  the  pole  In  the  s  plane,  then  far  from  the  surface  of  the  sphere 


the  modes  have  the  dependence 
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1 


nk 


(t  -  r/c) 


c  Z£k 


(42) 


nk 


Re(snk)-  ° 


where  we  are  using  the  traditional  electrical  engineering  convention 
and  s-*-*-jco,  and  the  complex  conjugate  appearing  in 
formula  (42)  is  used  to  agree  with  this  convention.  The  exponential 
form  of  (42)  has  the  well  known  "exponential  catastrophe"  behavior  as 
r  -*■  00  [22].  The  wave  amplitude  increases  in  the  radial  direction.  This, 
as  is  well  known,  is  accounted  for  in  setting  up  the  proper  excitation 
conditions  [7],  Creeping  waves  on  the  other  hand  obey  radiation 
conditions  at  infinity. 

An  interesting  question  that  arises  is  "how  do  the  pole  plots 
of  s  .  move  in  the  s  plane  when  the  object  geometry  changes  con- 

E1K 

tinuously?"  The  interpretation  of  the  resonance  conditions  (13)  and 
(15)  in  terms  of  the  residue  theorem  tells  us  that  the  pole  positions 
will  move  little  and  their  number  will  be  conserved  unless  an  addi¬ 
tional  singularity  of  the  field  is  created  by  the  geometry  deformation. 

A  deeper  understanding  of  this  analytic  function  method  should  be  pur¬ 
sued.  The  significance  of  topological  invariants  of  ("compact 
orientable  geometric")  surfaces  such  as  the  Euler-Poincare  character¬ 
istic  [18]  also  needs  to  be  investigated. 

An  immediate  logical  extension  of  this  theory  is  to  properly 
define  the  independent  closed  paths  f  for  more  general  surfaces. 

The  theory  should  then  be  checked  against  the  known  results  for  the 
prolate  and  oblate  spheroid. 
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APPENDIX  A 


The  solution  to  the  transcendental  equation 
3/2  . 


2/3  C 


Oz|!ll/2  dt 


(A.  1) 


is  now  given.  Olver  has  determined  chat  "where  branches  are  defined 
to  take  their  principal  values  when  z£(0,l)  and  ?£(0,°°)  and  are 
continuous  elsewhere"  (p.  421  of  [19]) • 

The  integral  is  doable: 


2/3? 


3/2  _ 


(A. 2) 


To  aid  in  the  solution  of  this  implicit  function  (i.e.,  given 
?  find  2)  let  us  make  the  change  of  variable 


z  *  secha  then  2/3?^  *  a  -tanh  a  (A.  3) 

Thus,  given  ui  where 

<u  *  a  -  tanh  a 
a  is  determined. 


Let  a)  ■  u  +  iv,  a  ■  a  +  i@  ,  tanh  cr  »  t  +  iy 

where  T  *  sinh2a/cosh(2a)  +  cos (2$)) 

U  -  sin  20/(cosh(2a)  +  cos(2g)) 


cosha  sing  -  i  sinh  a  cos8 
Notice  then  - - sinh‘a  +  sin‘S 


so  that 


u  ■  a 


sinh2ct 

cosh2a  +  cos 22 


(a) 


(A.  4) 


v  »  2 - - *5  (b) 

cosh2a  +  cos22 


Thus,  ia  (A. 4)  given  u  and  v,  ve  must  solve  for  (a, 2)  simultaneously. 


cosh2a  +  cos22  »  sinh2a  ^  cosh2a  +  cos22  *  ~ —ft 

Qe  ■*  U  P  —  V 

Therefore,  3inh2a  >  sin  28 

P“V 


Procedure : 

Given  (u,v)  cake  approximate  2  value  and  solve  for  correspond¬ 
ing  a  in  (A. 5).  Then,  take  this  (a, 2)  pair  and  substitute  them  into 
(A. 4a).  Vary  2  with  fixed  a  until  (A. 4a)  is  satisfied.  Take  this  new 
value  of  2  and  substitute  it  into  (A. 5)  to  obtain  an  updated  value  of 
a.  Stop  iteration  when  j  -  zn_^j  <  £  where  £  is  preset  tolerance. 


